Method for optimizing power production in photovoltaic modules

ABSTRACT

A method allows positioning the tracker at angles that promote the cooling of photovoltaic modules and therefore, decrease their operating temperature, without reducing the total energy produced thus optimizing power production of photovoltaic (PV) electricity by reducing the working temperature of PV modules of a solar tracker. The method also provides an optimization of the electrical output of the system for particular conditions of instantaneous air temperature and wind speed to improve electrical power generation ratios with respect to those known current techniques taking into account incident power in the PV plane, or in some cases the output power without considering the action of changes in wind speed or air temperature.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a national stage under 35 U.S.C. § 371 of PCT patent application PCT/ES2021/070197 filed on 18 Mar. 2021, which is hereby incorporated by reference in its entirety for all purposes. PCT/ES2021/070197 claims priority to European Patent Application 20169483.3 filed on 14 Apr. 2020, which is hereby incorporated by reference in its entirety for all purposes.

TECHNICAL FIELD

It is desirable to develop fast and effective methods that may optimize the production of electrical power without particularly requiring the operating temperature of the module, avoiding trial and error processes currently known in the art to determine the optimum operating point. Besides, this should come along with an accurate estimation of the electrical power that a photovoltaic (PV) module must produce under real operating conditions, improving the identification of failures and improper operational modes of the system by comparison between estimated power and produced power.

BACKGROUND

Currently, the module temperature variable is not considered when positioning solar panels, in fact it is assumed to be constant when the angle of orientation of the panel changes. The investigations carried out indicate that the cooling of the tracker, and therefore the operating temperature of the module, does vary according to the angle of position of the tracker. These variations have an influence on the production and, therefore, the consideration of this variable in the positioning of the tracker can increase the electrical production of the module.

Known methodologies that optimize electrical power production are mainly based on measurements without previously estimating the module's operating temperature and are based on trial and error processes to determine the optimum operating point. In addition, it avoids the indeterminacy situation under new conditions, in which the algorithms based on measurements take unknown positions. In addition, the complete simulation of the plant considering the topography allows reaching a combination of tracker positions that enables the optimization of global energy generation, instead of independent optimization of trackers.

DESCRIPTION

The object of invention allows to position the tracker at angles that promote the cooling of the modules and therefore, decrease their operating temperature, without reducing the total energy produced thus optimizing power production of photovoltaic (PV) electricity by reducing the working temperature of PV modules of a solar tracker. This, as the skilled person would acknowledge, increases the power conversion ratio enhancing the electrical conversion efficiency of the module in operation and, therefore, increasing electrical power production. The object of invention does not only solve the problems known in the art, but also provides an optimization of the electrical output of the system for particular conditions of instantaneous air temperature and wind speed thus improving electrical power generation ratios with respect to those known current techniques taking into account incident power in the PV plane, or in some cases the output power without considering the action of changes in wind speed or air temperature.

The method of the invention is preferably implemented in the tracker control system. The control system processes the input data and determines the angle at which the solar trackers should be positioned. In this way, the photovoltaic modules are oriented in the plane indicated by the algorithm, in which the electrical production is optimized, considering both the power received in the form of irradiation and the operating conditions involved in the power conversion of the module.

The method of the invention achieves a reduction in the temperature of the PV modules compared to traditional methods that simply arrange the module to maximize power production in the form of radiation received in the plane of the module. The method of the invention accounts for the fact that PV module temperature influences the efficiency of the photovoltaic conversion, and therefore, the higher the module temperature, the less electrical power produced.

Hence, the proposed method manages to estimate the operating temperature of the PV module as a function of the angle of movement of the tracker, and move the PV module to positions where, considering the power also received as radiation, the module temperature is lower, such that it improves power conversion efficiency and maximizes electrical production.

Although there are no similar methods, the method of the invention represents an improvement in reliability and speed compared to methods that use power measurements to optimize the power generated. The method of the invention determines and measures environmental variables such as air temperature and does not require access to power monitoring data produced by the plant. This represents an advance with respect to methods that optimize electrical production with all the intervening variables added, because it avoids a data acquisition system with sensors in the power cables, or feedback from the monitoring system outside the control system of the solar tracker. In this way, inputs of measured variables are reduced, and reliability of the positioning method is thereby increased and can be implemented by only considering air temperature value predictions.

The immediacy of the proposed method is effective from the first day, and therefore, avoids needing a learning period by overcoming this drawback of machine learning methods. The learning period of the machine learning methods for positioning the tracker that considers air temperature as the input variable is usually quite long, since the temperature presents seasonal variation throughout the year, which, moreover, is not directly coupled with irradiation. The proposed method is adjusted prior to the start of operation of the plant and does not require this learning period. The proposed method is also immune to initial malfunctions in the production of the solar plant. It is worth mentioning that in the case of PV plants, the modules usually have an initial adjustment period in which they may suffer abnormal initial degradations, which lengthen or introduce errors in the learning period of machine learning methods. Finally, the nature of the method, and the fact of having the air temperature decoupled, avoids the erroneous movement of the tracker under anomalous conditions that have not occurred previously, such as heat or cold waves, which are increasingly frequent because of climate change.

This method makes a prediction of the temperature of the solar module based on a physical model and moves the tracker directly to an optimal position; thus being faster and more effective than the observation and disturbance methods, avoiding losses associated with the movements of disturbance to different positions that may or may not be more optimal, as well as the observation period during which it is not in the sought position.

The analysis of the actual production data and the comparison with the estimated value may also be applied in the identification of failures and malfunction of the PV generation system. In this sense, as already described, the method of the invention achieves a reduction in the temperature at which the PV modules work, to increase the power conversion efficiency of the module in operation and therefore, an increase in the electrical power produced. It is much faster and more effective than methods that optimize electrical power without particularly considering the module operating temperature, since it avoids trial and error processes to determine the optimum operating point. It also allows a precise estimation of the electrical power that a PV module must produce under real operating conditions, so that it improves a process of identifying faults and malfunctioning of the system by comparing estimated power with produced power.

In order to do so, the object of the invention estimates the positioning angle of the PV modules at which electrical power produced is deemed to be at its maximum value which depends on the input power to the panel in the form of radiation, and on the power conversion efficiency of the panel, which in turn depends on the panel operating temperature. The proposed method estimates the module's operating temperature as a function of the module's angle, air temperature, wind speed, the albedo of the soil, the type of panel and the geometric conditions of the distribution of the trackers (distance between corridors, layout, topography). For each position of the sun, the method is applied to the conditions of each solar tracker and estimates the power generated considering both the power received as a function of the angle and the operating temperature and decides which is the most optimal angle.

A complete consideration of the solar plant including the topography reaching a combination of tracker positions that yields a maximized global power generation, instead of an independent optimization of trackers. This aspect is hardly achievable by methods based on measured data or machine learning, since it would be unfeasible to consider the range of possible combinations that the application of predictive methods allows.

For this, the method of the invention implements a loss power estimation function based on the angle of inclination that balances with the temperature, obtaining an optimal angle. Additionally, a time buffer that analyzes the variability based on the wind speed to avoid continuous or alternated movements of the tracker may be implemented.

BRIEF DESCRIPTION OF FIGURES

To complement the description that is being made and to aid towards a better understanding of the characteristics of the present invention, in accordance with a preferred example of practical embodiment thereof, a set of drawings is attached as an integral part of said description wherein, with illustrative and non-limiting character, the following has been represented:

FIG. 1 shows a diagram showing the steps of the method of the invention

FIGS. 2 a, 2 b show respective diagrams where the definition of the different angles mentioned hereby angles are represented; therefore, providing an overview of a tracker orientation angle modification procedure.

FIGS. 3 a, 3 b respectively depict diagrams showing a general description of the tracker orientation angle modification method assessing power losses due to misorientation for a constant temperature T=25°.

PREFERRED EMBODIMENT

The advantages of the invention are better understood from FIG. 1 where the object of the invention is illustrated as a flowchart diagram.

As an added advantage of the method, it allows a precise estimation of the electrical power that a PV module must produce in real operating conditions, so that it improves a process of identifying faults and malfunctioning of the system by comparing estimated power and power produced.

As per FIGS. 2 a and 2 b , the method of the invention encompasses an estimation of the module operating temperature which is made for a battery of possible tracker angles, based on a set of predetermined variables. These temperatures values are applied to the power received in the form of radiation (Eα) applying linearly the temperature coefficient at maximum power (γ) defined by the panel manufacturers, determining the Efficient power or power affected by the variation in the efficiency of conversion (Ef). Finally, the angle that maximizes this Efficient power is selected as the angle for the orientation of the tracker.

The module temperature for each angle of the tracker is calculated using a thermodynamic model of the tracker based on environmental variables to which the module is exposed, which incorporates model adjustment parameters that describe the tracker and plant geometry, thus such as the thermodynamic and power properties of the module.

-   -   The environmental input variables of the model may account as:         Air temperature, wind speed W and angle of wind direction σ,         beam irradiation Gbeam and position of the sun, diffuse         irradiation Gdiff and the albedo of the soil Albedo. Note that         beam irradiation Gbeam or Direct Normal Irradiance (DNI) is         measured at the surface of the Earth at a given location with a         surface element perpendicular to the Sun excluding diffuse solar         irradiation (irradiation that is scattered or reflected by         atmospheric components). Direct irradiance is equal to the         extraterrestrial irradiance above the atmosphere minus the         atmospheric losses due to absorption and scattering and losses         depend on time of day (length of light's path through the         atmosphere depending on the solar elevation angle), cloud cover,         moisture content and other contents. The irradiance above the         atmosphere also varies with time of year (because the distance         to the sun varies), although this effect is generally less         significant compared to the effect of losses on DNI. Diffuse         irradiation Gdiff, also called sky irradiation, diffuse         skylight, or just skylight, is solar irradiation reaching the         Earth's surface after having been scattered from the direct         solar beam by molecules or particulates in the atmosphere.     -   The geometric parameters of the plant and the tracker may         account as: Distance between trackers or pitch, topography,         tracker height, separation between panels, dimensions and gaps         in the tracker structure.     -   The parameters derived from the properties of the module, the         nominal power conversion efficiency of the module, the         bifaciality (for monofacials this value is close to zero), color         and characteristics of the surface, type and shape of the frame.

The method produces results based on a balance between:

-   -   Power falling upon the panel as solar irradiation (front and         rear) being absorbed by the module depending on the properties         of its surface. This part includes the calculation of the         irradiation for the surface of each tracker considering the         topography and North-South inclination of each tracker.     -   Power delivered in the form of electricity produced (power         conversion efficiency) by both front and rear parts, namely         power efficiency allowing determining a percentage of         irradiation converted into electricity.     -   Power transferred to the air and to the rest of the bodies in         the scene as heat (refrigeration) namely, an amount of heat         transferred by refrigeration is estimated considering         convective, conductive and radiation factors also accounting for         geometric properties of the plant and the topography as well as         the position of the set of trackers in the plant. The effect of         positions of adjacent trackers which may influence temperature         values of other related trackers may be considered as well.

The method may also incorporate a time buffer that analyzes the variability based on the historical wind speed to obtain a stable parameter thus avoiding contiguous position corrections due to gusts.

The complete simulation of the plant considering the topography allows reaching a combination of tracker positions that enables to maximize global generation, instead of independent optimization of trackers. This aspect is hardly achievable by methods based on measured data or machine learning, since it would be unfeasible to consider the range of possible combinations that the application of predictive methods permits.

The main advantage is that it optimizes the electrical production of the system for conditions of instantaneous air temperature and wind speed. It improves the electrical generation with respect to the current technique, considering incident power in the PV plane, or in some cases the output power without considering the action of changes in wind speed or air temperature.

The goal of the method of the invention is based on the balance of the losses due to temperature and/or misorientation. In a preferred embodiment of the object of the invention depicted in FIGS. 3 a and 3 b , power due to misorientation for a constant module temperature T=25° is assessed: Pα=E(β)−E(α) wherein: E(α)=μ(Gdiff(α)+Gbeam(β)cos(β−α)+bifaciality·Grear(α)) E(β)=μ(POA+bifaciality·Grear(β))

being:

Pα=power penalty function due to misorientation

t=time and date

Eα=module electrical power estimation at a given temperature value, in this preferred embodiment T=25°

POA=Plane sun-oriented power being the power in the form of irradiation that reaches the front plane of the module.

α=Non sun-oriented tracker angle

β=Optimal sun-oriented tracker angle

Bifaciality=module bifaciality (0≅1)

Albedo=soil albedo

μ=PV system electrical conversion efficiency

Grear=Backplane irradiation

Gbeam=Front plane beam irradiation

Gdiff=Front plane diffuse irradiation

The object of the invention is preferably based on on-field measured irradiation values, said values may be collected by sensors configured to capture, among others, backplane irradiation (Grear), front plane beam irradiation or (Gbeam) or Front plane diffuse irradiation (Gdiff) values.

In some embodiments of the invention module temperature as a function of the angle T(α) may be assessed as:

T(α) = (k₀ + k₁T_(air) + +(k₂ + k₃p_(air)T_(air) + k₄W_(air)+ k₅(cos (σ)W + k₆cos²(σ)W²) ⋅ (k₇ + k₈d + k₉d²) + + $\sum_{i = 0}^{4}{{{k_{10 + i}\left( {{{albedo} \cdot {POA} \cdot d \cdot \cos}\alpha} \right)}^{i}++}{k_{11}\left( {1 - \mu} \right)}\left( {{{POA}\left( {D + 1 - \frac{\left( {\alpha - \beta} \right)^{2}}{2}} \right)} +} \right.}$ Grear(k₁₂ − k₁₃bfacilaity)) + k₁₄T(α)^(1/4))

Being:

p_(air) atmospheric pressure

T_(air) Air temperature

d pitch or distance between trackers

W Wind speed

σ angle of wind direction

μ solar module efficiency

θ angle of the solar tracker

k_(i) adjustment parameters of the solar tracker thermodynamic model.

D=Relationship between beam and diffuse irradiation

Wherein:

T_(air) is preferably measured in the field by sensors, eventually it may also be determined based on weather forecast or satellite imaging-based predictions.

Wind speed and direction are preferably measured on site by means of one or more wind vane anemometers.

In yet an alternative embodiment, another penalty function may be implemented accounting for the power effort of moving the tracker and the weight of inconvenience of investing time and how far is the extreme point of minimum, as well as an estimation of the probability that the wind changes the next minutes and the future temperature. This penalty function is used in order to decide if it is interesting to carry out the change. Pm(α)=m ₁(θ−α)w ²

being:

Pch=Penalty function evaluating the power losses due to angle change.

αo=current Tracker angle

Prt(α, αo)=power used for rotating the tracker

Pmv(α, αo)=power losses during the tracker movement due to pass lower irradiation points

Pwp(α, αo)=power losses temperature in case of future wind speed with a probability value around 20%.

Being θ as the positioning angle of the tracker θ=β, the method of the invention comprises:

Calculation of an interval of possible a in a range between −15° and +15° around the Oriented angle doing so in one-degree steps. α=β−Δ . . . β+Δ Calculation of disorientation penalty when going to an angle P (α), as per: P(α)=E(β)−E(α) Calculation of the temperature gain Pt (α), based on the updated environmental variables as per: Pt(α)=E(β)γ(T(β)−T(α)))

Obtaining a value for angle α that minimizes the loss due to disorientation Pα−Pt(α).

Evaluation of the loss due to disorientation for said angle that minimizes the loss due to disorientation:

-   -   Should the value of said loss determined to be positive, the         tracker is actuated to go to the oriented position, hence θ=β.     -   Should the value of said loss determined to be negative, the         method comprises evaluating penalty functions:         Pm(α)=m ₁(θ−α)w ²         Pch(α)=c ₁ w ² +c ₂ cos cos σ         P(θ)=E(β)−E(θ)         Pt(θ)=E(β)γ(T(β)−T(θ)))

Then, a calculation of penalty for movement between current position and the estimated Pm (α), and penalty for confidence in the measure Pch (α) disorientation between the current position and the oriented P (θ) and losses due to temperature gain between the current position and the oriented Pt (θ) are carried out.

Evaluation of the sum of the above terms together with the disorientation penalty P(α) and temperature gain Pt (α) from the angle β as per: +P(α)−P(θ)+Pt(α)−Pt(θ)+Pm+Pch

If this resulting value is positive, the tracker is not actuated, therefore θ=θ

If this resulting value is negative, the tracker is actuated moving to the estimated position where the temperature optimum angle is achieved, therefore θ=α and the method starts over again from the calculation of an interval of possible α in a range between −15° and +15° around the Optimal sun-oriented tracker angle (β) preferably doing so in one degree steps. α=β−Δ . . . β+Δ 

The invention claimed is:
 1. Method for optimizing power production in photovoltaic modules wherein θ is an angle of the tracker, the method comprising determining a value of the following variables: a. air temperature (T_(air)), b. atmospheric pressure (p_(air)), c. wind speed (W), d. angle of wind direction (σ), e. plane sun-oriented power (POA) being the power in the form of irradiation that reaches the front plane of the module, f. backplane irradiation (Grear), and g. relationship (D) between Gbeam=Front plane beam irradiation, and Gdiff=Front plane diffuse irradiation, the method being characterized by comprising: i. calculation of an interval of possible values for a non sun-oriented angle (α) comprised in a range between −15° and +15° around an optimal sun-oriented tracker angle (β) wherein irradiation value maximizes, ii. calculation of disorientation penalty when moving the tracker to an angle P (α), as per: P(α)=E(β)−E(α) wherein: E(α)=μ(Gdiff(α)+Gbeam(β) cos(β−α)+bifaciality·Grear(α)) E(β)=μ(POA+bifaciality·Grear(β)) Pα=power penalty function due to misorientation, Eα=module electrical power estimation at a given temperature value T, POA=plane sun-oriented power being the power in the form of irradiation that reaches the front plane of the module, α=non sun-oriented tracker angle, β=optimal sun-oriented tracker angle, assessed using any of method from the state of the art, in a preferred embodiment of the invention astronomic tracking techniques were used, Bifaciality=module bifaciality μ=PV system electrical conversion efficiency, Grear=backplane irradiation, Gbeam=Front plane beam irradiation, and Gdiff=Front plane diffuse irradiation calculation of the temperature gain Pt (α), based on updated environmental variables selected from the group consisting of: Air temperature (T_(air)), wind speed (W) and angle of wind direction (σ), beam irradiation and position of the sun, diffuse irradiation and albedo of the soil; as per: Pt(α)=E(β)γ(T(β)−T(α))) wherein: T(α) = (k₀ + k₁T_(air) + +(k₂ + k₃p_(air)T_(air) + k₄W_(air)+ k₅(cos (σ)W + k₆cos²(σ)W²) ⋅ (k₇ + k₈d + k₉d²) + + $\sum_{i = 0}^{4}{{k_{10 + i}\left( {{{albedo} \cdot {POA} \cdot d \cdot \cos}\alpha} \right)}{i++}{k_{11}\left( {1 - \mu} \right)}\left( {{{POA}\left( {D + 1 - \frac{\left( {\alpha - \beta} \right)^{2}}{2}} \right)} +} \right.}$ Grear(k₁₂ − k₁₃bifaciality)) + k₁₄T(α)^(1/4)) being: p_(air) atmospheric pressure T_(air) Air temperature d pitch or distance between trackers W Wind speed σ angle of wind direction μ solar module efficiency θ tracker angle k_(i) adjustment parameters of the tracker thermodynamic model iii. obtaining a value for the non sun-oriented angle (α) that minimizes losses due to disorientation Pα−Pt (α), iv. evaluation of the loss due to disorientation for said non sun-oriented angle (α) that minimizes losses due to disorientation Pα−Pt (α) in such a way that: a. should the value losses due to disorientation Pα−Pt (α) determined to be positive, the tracker is actuated to go to the oriented position related to an angle wherein front plane irradiation level Gbeam is maximum, hence θ=β b. should the value of said losses due to disorientation Pα−Pt (α) determined to be negative, the method comprises evaluating: penalty for movement Pm(α) between current position and the estimated position following: Pm(α)=m ₁(θ−α)w ², being m₁ a parameter of the thermodynamic model, penalty for confidence in measurements Pch(α) following: Pch(α)=c ₁ w ² +c ₂ cos cos σ being c₁, c₂ parameters of the CFD model, disorientation penalty P(θ) between a current position and an oriented position following: P(θ)=E(β)−E(θ) temperature gain Pt(θ) between the current position and the oriented position: Pt(θ)=E(β)γ(T(β)−T(θ))) c. evaluation of the losses together with the disorientation penalty P(α) and temperature gain Pt (α) from the optimal sun-oriented tracker angle (β) as per: +P(α)−P(θ)+Pt(α)−Pt(θ)+Pm(α)+Pch(α) in such a way that: resulting value being positive, the tracker is not actuated, therefore θ=θ, and resulting value being negative, the tracker is actuated moving to the estimated position where a maximum energy production level considering temperature optimum angle is achieved, therefore θ=α; and the method comprises starting over again from step i.
 2. The method of claim 1 wherein the air temperature (T_(air)) is determined by means of sensors associated to the panel.
 3. The method of claim 1 wherein the wind speed (W) is determined by means of sensors associated to the panel.
 4. The method of claim 1 wherein the angle of wind direction (σ) is determined by means of sensors associated to the panel.
 5. The method of claim 1 wherein at least one irradiation value selected from the group consisting of: backplane irradiation (Grear), front plane beam irradiation or (Gbeam) or Front plane diffuse irradiation (Gdiff), is picked up by on-field sensors. 